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ABSTRACT  AND  INTRODUCTION 

Since  1978,  there  has  been  a  history  of  successful  cooperation  between  control 
researchers  at  the  University  of  Minnesota  and  the  Hungarian  Academy  of  Sciences. 
This  cooperative  research  effort  was  initially  sponsored  by  the  National  Science 
Foundation  in  the  late  1970's  and  early  1980's  and  later  it  has  continued  on  an  informal 
basis  to  the  present. 

One  of  the  most  important  areas  of  research  in  control  theory  is  the  design  of 
feedback  controllers  for  systems  which  have  significant  uncertainties  in  the  plant.  These 
uncertainties  can  result  from  a  lack  of  precision  in  mathematical  modeling  of  the  plant 
and/or  changes  in  the  plant  parameters  with  time.  Two  main  techniques  for  design  of 
controllers  for  systems  with  significant  uncertainties  are  robust  and  adaptive  control 
theories.  System  identification  is  an  important  tool  in  both  techniques  as  well  as  in 
many  other  control  design  methods. 

The  US  cooperating  partner  has  a  long  experience  in  designing  high  performance 
robust  control  systems  (helicopter  flight  control,  supermaneuverable  aircraft  control, 
missile  autopilots,  etc.).  The  Hungarian  partner  has  implemented  several  high 
performance  adaptive  controllers  for  complex  processes  (national  load-frequency 
control  system,  chemical  composition  control  at  cement  raw  material  handling, 
combined  fineness  and  effectiveness  control  at  closed  circuit  grinding  mills,  etc.).  Both 
sides  have  long  experience  in  modeling,  structure  and  parameter  estimation  of 
sophisticated  dynamic  systems  and  have  obtained  significant  achievements  in  robust 
control  methodologies.  Therefore  a  fruitful  cooperation  was  expected  from  a  new 
planned  joint  research  effort,  targeting  the  complex  problems  of  simultaneous 
identification  and  control,  furthermore  on  the  interaction  between  systems  identification 
and  robust  and  adaptive  control  theories.  Applications  could  include  helicopter  flight 
control  and  active  vibration  control. 

During  the  preparation  of  the  project  and  the  evaluation  of  the  preproposals  the 
European  Research  Office  of  the  US  Army  suggested  that  it  was  desirable  to  submit  a 
separate  proposal  by  the  Hungarian  side  under  the  title  "Identification  For  Robust 
Control",  because  the  ARO  is  not  prepared  yet  to  handle  the  above  type  of  joint 
projects.  Finally  this  separate  proposal  was  accepted  requesting  yearly  evaluation.  So 
this  is  the  final  report  of  the  first  year's  activity  of  the  Hungarian  side.  Besides  the  US 
and  Hungarian  partners  promised  to  keep  each  other  informed  about  their  achievements 
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and  try  to  form  a  loosely  formalized  connection  according  to  the  original  ideas. 

KEYWORDS 

System  identification;  process  parameter  estimation;  optimal,  robust  and  adaptive 
filtering  and  control;  combined  identification  and  control;  identification  for  control. 
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PROJECT  SUMMARY 

The  research  contract  was  signed  on  March  24,  1994,  so  the  reported  period  is  about  15 
months.  This  period  was  enough  to  achieve  significant  and  completed  contributions.  In 
this  reported  period  we  could  make  a  very  intensive  research  activity  in  the  areas 
indicated  and  planned  in  the  project  proposal  shown  by  the  large  number  of  papers. 

Herebelow  we  try  to  summarize  only  those  achievements  that  we  could  formalize  in 
written  publications  or  reports.  Because  most  of  them  are  submitted  papers  to 
international  scientific  journals  and  conferences,  there  is,  of  course,  a  certain  chance 
that  not  all  of  our  materials  will  be  accepted.  However,  we  hope  that  most  of  them  will 
finally  be  published,  because  our  group  has  a  good  acceptance  rate. 

Furthermore  our  recently  published  and  submitted  papers,  initiated  and  supported  by 
the  project,  are  also  summarized  to  show  the  development  of  our  work.  Those  working 
papers  (project  reports)  are  also  presented,  which  are  for  preparing  and  founding  our 
further  research. 

We  planned  our  research  activities  basically  in  three  directions  : 

(1)  the  general  (classical?)  robust  control  approach 

(2)  special  auxiliary  solutions  augmenting  robustness  of  adaptive  control 

(3)  a  new  approach  based  on  a  generic  optimal  control  system  structure. 


The  general  robust  control  approach 

In  this  area  we  wanted  to  follow  the  developments  in  the  classical  robust  control 
approaches  and  to  reach  our  new  contributions  along  these  lines.  (The  technical  details 
of  our  expectations  what  and  how  to  do  in  our  planned  research  were  given  in  our 
project  proposal.) 

H2  and  identification  for  control 

We  introduced  a  new  interpretation  of  criteria  in  frequency  domain.  This 
interpretation  gives  special  "butterfly"  sensitivity  (or  robustness)  shaping  regions  that 
can  easily  be  applied  on  the  classical  Bode  plot  design  techniques  [9],  [19]  and  [27]. 
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This  interpretation  seems  to  be  very  useful  to  handle  the  critical  medium  frequency 
domain  in  combined  or  simultaneous  identification  and  control  problems. 

Various  approaches  to  identification  of  dynamic  systems  using  criteria  have 
been  analyzed.  The  report  [16]  provides  a  comprehensive  survey  and  evaluations  of  the 
linear  and  nonlinear  methods  based  on  interpolatory  algorithms.  Numeric  realizations 
of  these  methods  are  also  provided  in  MATLAB™  environment  and  this  can  serve  to 
make  comparisons  using  either  simulated  or  measured  data. 

In  [20]  and  [29]  we  considered  the  identification  of  systems  with  transfer  function  in 
the  disc  algebra.  The  approximate  model  is  a  finite  n  -order  FIR  model.  Using  linear 
least  squares  parameter  estimation  with  a  properly  designed  input  signal,  it  is  shown 
that  the  asymptotic  estimates  obtained  from  time  domain  data  are  the  first  n  Taylor 
coefficients.  This  model  has  the  minimum  mismatch  error  norm  over  all  FIR 
models  that  are  polynomial  in  the  delay  operator  with  order  not  exceeding  n . 

There  is  a  renewed  interest  in  the  choice  and  selection  of  appropriate  basis  for 
identification.  This  is  illustrated  by  the  fact  that  special  sections  have  been  organized  in 
the  SYSID’94  and  CDC'94  symposia.  In  [12]  we  showed  that  a  basis  in  the  disc  algebra 
can  be  obtained  from  the  Faber-Schauder  and  Franklin  systems.  The  dynamical  systems 
with  transfer  functions  in  the  disc  algebra  can  be  approximated  in  norm  by  partial 
sum  operators.  We  elaborated  a  linear  algorithm  for  computation  of  approximate 
models  from  frequency  domain  data.  Numeric  computations  and  comparisons  to 
existing,  like  Cesaro-mean  based  linear  algorithms  showed  a  faster  convergence  with  a 
smaller  error  norm. 

Convergence  analysis  for  deterministic  case  and  under  "unknown  but  bounded" 
noise  assumptions  are  provided  in  [17]. 

and  robust  filtering 

The  design  of  filters  with  performance  appears  in  the  output  estimation  problem 
in  relation  to  the  output  feedback  control.  The  motivation  is  to  design  a  filter  or 
observer  that  guarantees  a  given  attenuation  of  an  external  disturbance  signal  that  also 
appears  in  the  filter  state  or  output  error  system. 

The  robustness  problem  is  related  to  the  situation  where  the  filter  is  designed  on  the 
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basis  of  a  nominal  system  while  there  is  an  additive  or  multiplicative  error  transfer 
function  between  this  nominal  transfer  function  and  the  transfer  function  of  the  real 
system. 

The  usual  solution  to  this  robustness  problem  that  has  been  published  in  the  past  few 
years  has  a  dual  structure  when  compared  to  the  two-Riccati  equation  type  state  space 
solution  of  the  output  feedback  control. 

In  our  approach  we  considered  dynamic  systems  given  by  a  LTI  nominal  model  and 
by  time-varying  parametric  perturbations  in  the  state  space  matrices. 

Using  a  basis  of  the  Lie-algebra  generated  by  the  nominal  state  transition  matrix  and 
parametric  disturbance  matrices  we  construct  an  "auxiliary"  LTI  system  that  has 
invariant  subspaces  related  to  parametric  perturbations  that  are  relevant  in  filter  design 
that  are  identical  to  the  invariant  subspaces  of  the  original  LTV  system.  The  effect  of 
parametric  uncertainties  appears  in  this  auxiliary  system  as  unknown  but  bounded  £2 
externals  signals.  Using  this  approach  the  problem  can  be  reduced  to  the  design  of  a 
filter  to  attenuate  the  effect  of  the  external  signal  in  the  state  or  output  error  in  L2  norm 
as  much  as  possible  (worst  case  disturbance  concept),  leading  to  the  design  of 
filters. 

This  approach  was  applied  to  the  design  of  fault  detection  filters  where  the  alarm 
levels  of  the  filters  can  be  adjusted  to  the  maximum  attenuation  of  system  uncertainties 
represented  by  external  deterministic  noise  in  the  auxiliary  system.  Results  associated 
to  this  topics  can  be  found  in  [21], [32], [34]  and  [36]. 

Our  further  contributions  are  still  remarkable  in  designing  robust  diagnostic  and 
failure  detection  systems  that  are  based  on  special  classes  of  system  identification 
methodologies  [6], [7], [15]. 


Special  auxiliary  solutions  augmenting  robustness  of  adaptive  controls 

In  this  area  we  investigated  special  methods  that  augment  the  robustness  of  classical 
(mostly  adaptive)  control  systems.  Our  contributions  are  remarkable  on  the  area  of 
adaptive  time  delay  estimation  [1],[2],[4].  These  methods  considerably  extend  the 
applicability  of  adaptive  regulator  algorithms. 
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Besides  investigating  adaptive  time  delay  estimation  eonnected  to  PID  regulators, 
the  robustness  of  simple  PID  controls  has  also  been  treated.  A  new  design  principle  was 
introduced  decreasing  the  gain  sensitivity  of  these  regulators  [11], [25].  New  adaptive 
PID  regulator  algorithms  were  also  developed  as  virtual  instruments  (National 
Instruments'  LabView™  environment)  to  reach  easy  portability  and  wide  applicability 
[3],[5]. 

On  the  basis  of  these  theoretical  results  and  simulation  studies  we  start  to  investigate 
their  applicabilities  for  helicopter  models  targeting  to  develop  appropriate  adaptive  PID 
regulators  for  their  control.  These  applications  can  also  be  interesting  for  ARO. 


A  new  approach  based  on  a  generic  optimal  control  system  structure 

We  believe  that  one  of  the  major  theoretical  contribution  of  our  project  will  be  a  new 
regulator  design  approach  based  on  a  generic  optimal  control  scheme  developed  by  us. 
This  is  a  new  structure  how  to  design  optimal  pole  placement  controllers.  This  new 
scheme  allows  to  avoid  the  explicit  solution  of  a  polynomial  equation  obtaining  the 
transfer  function  elements  of  the  optimal  controller  directly  [8].  Our  recent  papers  [18], 
[22],  [23],  [28]  and  [31]  discuss  the  different  variants  and  details  of  this  new  design 
methodology.  Many  further  publications  are  planned  for  the  different  applications  of 
these  principles. 

The  new  design  principle  is  quite  general  and  applicable  for  nonminimum  phase 
(inverse  unstable)  and  delay  time  systems,  too.  The  controller  is  easy  to  be 
implemented  in  computer  controlled  systems. 

The  structure  of  the  optimal  controller  gives  a  special  insight  to  understand  the 
operation  of  a  feedback  loop  for  the  servo  and  disturbance  rejection  paradigm.  It  is  easy 
to  see  the  role  of  the  knowledge  of  the  system  model,  the  role  and  appearance  of  those 
system  factors  that  are  invariant  to  any  control  strategy. 

This  generic  structure  opens  a  new  way  to  handle  system  uncertainties,  furthermore 
combined  modeling  and  control  issues  providing  a  new  canonical  sensitivity  scheme. 

We  formulated  a  special  framework,  where  the  different  separate  phase  of 
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identification  for  and  design  of  robust  control  can  properly  be  handled.  This  approach 
combines  the  classical  "minimum  variance"  like  control  with  the  new  concept  of 
"maximum  variance"  input  design  for  robust  identification  for  control.  New 
publications  are  under  preparation. 

This  generic  controller  structure  opens  a  new  way  how  to  handle  system 
uncertainties,  furthermore  combined  modeling  and  control  issues  providing  new 
canonical  sensitivity  schemes.  We  investigated  these  new  schemes  and  surprising  new 
results  have  been  obtained  [24]  and  [30]. 

The  properties  of  the  simple  combination  of  the  open-loop  identification  and  closed- 
loop  control,  the  combined  identification  and  control  in  a  parallel  closed-loop  optimal 
controller  scheme,  in  a  parallel  in-loop  optimal  controller  scheme  and  in  a  new  generic 
optimal  controller  scheme  were  investigated.  It  was  shown  that  the  new  generic  optimal 
controller  structure  is  superior  to  the  others. 

The  simple  combination  of  the  open-loop  identification  and  closed-loop  control  does 
not  provide  good  properties  for  iterative  control  refinement. 

The  widely  used  and  analyzed  parallel  closed-loop  optimal  controller  scheme  is 
good  for  iterative  control  refinement,  however,  it  really  needs  a  cautious  windsurfer 
approach,  because  its  robustness  is  highest  around  the  crossover  frequency  co^. ,  which 

is  not  known  apriori  and  can  be  approached  only  in  an  iterative  procedure. 

The  parallel  in-loop  optimal  controller  scheme  mostly  applied  at  adaptive  regulator 
schemes  is  directly  not  good  for  iterative  control  refinement,  having  improper 
frequency  domain  robustness  shaping.  It  is  however  possible  to  use  a  special  prefilter 
based  on  the  inverse  of  the  model  based  characteristic  equation  for  the  identification 
data  to  make  this  method  similar  to  the  previous  one. 

It  is  also  interesting  to  know  that  in  a  noise-free  case  the  last  two  methods  can  easily 
fail  because  the  identification  errors  can  be  zero,  if  the  control  errors  are  zero,  inspite  of 
a  nonzero  modeling  error.  This  bad  feature  can  stop  the  identification  at  a  not 
predictable  wrong  model  and  not  continuing  any  iterative  refinement. 

The  new  generic  optimal  controller  scheme  seems  to  be  best  from  the  investigated 
methods,  because  its  frequency  domain  robustness  shaping  amplifies  the  high 
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frequency  uncertainties  in  the  identification  error  and  more  accurate  model  is  expected 
in  the  targeted  high  frequency  domain  (where  we  want  to  increase  the  bandwidth).  It 
behaves  also  well  because  the  identification  and  control  errors  are  the  same  and  they 
can  vanish  only  at  the  same  time. 

In  the  first  interim  report  it  was  mentioned  that  we  formulated  a  special  framework, 
where  the  different  separate  phases  of  identification  for  and  design  of  robust  control  can 
properly  be  handled.  This  approach  combines  the  classical  "minimum  variance"  like 
control  with  the  new  concept  of  "maximum  variance"  input  design  at  robust 
identification  for  control.  In  that  phase  this  idea  was  only  a  conjecture.  Using  the 
CSILLA™  package  (below)  we  could  prove  the  nice  converging  properties  of  the  new 
"maximum  variance"  reference  input  design  algorithm  [10].  This  method  works  well 
both  in  the  iterative  off-line  and  adaptive  on-line  situation,  too.  The  very  interesting  and 
effective  feature  of  this  method  is  that  it  gradually  (iteratively  or  recursively,  depending 
on  the  applied  scheme)  improves  the  frequency  spectrum  of  an  initial  reference  input 
signal  excitation  approaching  and  concentrating  on  the  vital  medium  frequency  domain 
around  the  cross-over  frequency.  This  seems  to  be  a  natural  way  to  solve  the  well 
known  problem  of  "catch  of  22"  in  combined  identification  and  control  paradigms, 
when  one  should  know  the  process  to  generate  the  optimal  excitation  or  to  have  apriori 
or  "god  given"  frequency  dependence  of  modeling  error  sensitivities. 

This  new  principle  for  optimal  reference  input  design  outlines  a  new  scheme  of 
"trial-control",  which  is  a  special  possible  extension  of  the  classical  "dual-control" 
scheme  -  based  on  the  simultaneous  identification  and  control  loops  -  with  a  third  loop 
for  optimizing  the  excitation  [10],  [33].  Several  different  combinations  of  these  loops 
can  be  imagined  worth  studying  in  the  near  future. 

The  above  new  results  were  deeply  analyzed  and  compared  in  the  invited  survey 
paper  [10]. 


Elaboration  of  software  tools 

It  turned  out  that  we  need  special  software  tools  to  perform  "mathematical 
experiments",  i.e.,  special  simulation  runs  to  demonstrate  the  operation,  convergence, 
applicability  and  performance  behaviors  of  the  developed  methods  and  algorithms.  A 
careful  analysis  showed  that  the  MATIAB™  program  and  environment  widely  used  in 
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the  control  community  is  basically  good  for  our  goals,  if  we  write  a  couple  of  new 
routines  and  functions.  Another  point  was  that  in  this  way  we  are  compatible  with  our 
joint  research  partner  (U  of  M,  USA)  and  other  scientific  schools  in  this  area.  We 
finished  an  important  package  written  in  MATLAB™,  capable  to  handle  all  special 
cases  we  discussed  in  the  theoretical  papers  and  reports.  This  definitely  speeded  up  the 
experimental  proof  of  our  results. 

The  following  subsystems  have  been  elaborated: 

•  combined  identification  and  control  for  inverse  stable  processes 

•  combined  identification  and  control  for  inverse  unstable  processes 

•  combined  identification  and  control  for  inverse  unstable  processes  with 
optimal  reference  input  design 

•  adaptive  identification  and  control  for  inverse  unstable  processes 

•  adaptive  identification  and  control  for  inverse  unstable  processes  with 
optimal  reference  input  design 

•  comparison  of  step,  square-wave  and  the  new  maximum-variance  optimal 
reference  input  design  approach 

•  design  of  robust  BID  regulators 

•  combined  identification  and  design  of  robust  BID  regulators  with  optimal 
reference  input  design  based  on  the  new  "catamaran"  approach 

•  combined  identification  and  control  for  first  and  second  order  roll  response 
helicopter  model 

•  combined  identification  and  control  for  first  and  second  order  roll  response 
helicopter  model  with  optimal  reference  input  design 

•  adaptive  identification  and  control  for  first  and  second  order  roll  response 
helicopter  model 

•  adaptive  identification  and  control  for  first  and  second  order  roll  response 
helicopter  model  with  optimal  reference  input  design 

•  BID  control  for  first  and  second  order  roll  response  helicopter  model 

•  robust  BID  control  for  first  and  second  order  roll  response  helicopter  model 

The  above  newly  developed  package  of  subsystems  are  referred  as  CSILLA™ 
(Combined  System  Identification  and  controL  package  based  on  MATLAB™)  in  the 
following.  We  had  to  develop  special  gateways  to  the  well  known  general  word¬ 
processing  and  graphing  programs,  too,  to  ensure  publishing  quality  for  our  plots  and 
drawings.  We  used  the  Apple  Macintosh™  platform,  where  the  integrity  of  these 
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programs  is  provided  on  a  very  high  level.  Unfortunately  we  could  use  only  the  Mac  II 
version  of  MATLAB™  available  and  providing  an  acceptable  speed  for  our  interactive 
and  iterative  simulations.  However,  this  speed  is  still  very  far  from  the  originally 
required  high  performance  computing  facilities.  We  try  to  enhance  our  computing 
infrastructure  to  reach  the  PowerPCA^  or  workstation  performance  in  the  near  future. 

One  can  see  that  this  list  is  a  thorough  package  of  the  covered  areas.  We  mainly 
concentrated  on  the  developed  new  efficient  methods  for  combined  identification  and 
control.  These  are  mainly  off-line  methods  for  iterative  techniques.  The  adaptive 
versions  of  the  above  algorithms  have  also  been  elaborated. 


Application  possibilities 

We  started  to  study  vehicle  active  suspensions  based  on  H2/H^  optimization  robust 
control  theory,  which  can  have  special  interest  for  ARMY  applications  [13], [14]. 

We  started  to  investigate  the  applicability  of  these  new  theoretical  results  for 
improving  the  performance  of  possible  control  schemes  based  on  helicopter  models. 
We  used  theoretical  and  experimental  models  published  for  roll  response  modeling.  The 
roll  response  in  the  frequency  range  of  interest  is  essentially  second-order,  owing  to  the 
coupling  of  the  rotor  regressive  flapping  and  fuselage  roll  modes.  A  more  general 
second  order  model  involving  the  rotor  flap  inverse  time-constant,  the  total  flapping 
stiffness,  the  lateral  stick  gearing  and  the  residual  time  delay,  furthermore  a  simpler  first 
order  model  represented  by  a  time  delay  and  the  coupled  fuselage/rotor  dynamics  given 
by  the  roll  gain  and  the  first  order  roll  damping  were  considered. 

The  simulation  studies  using  these  helicopter  models  presented  how  to  increase  the 
robustness  of  PID  regulators  in  this  application  area.  The  new  advanced  methods  for 
combined  identification  and  control  of  these  models  were  also  treated.  Several  adaptive 
schemes  have  also  been  investigated.  It  was  also  shown  what  is  the  influence  and 
benefit  if  we  use  optimal  reference  signal  design  to  improve  the  convergence  of  the 
combined  and  adaptive  schemes. 

These  examples  were  used  as  simulation  examples  in  our  recent  publications, 
however,  we  plan  to  write  a  special  report  on  these  application  possibilities  and  results 
in  the  second  year  of  the  project. 
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During  the  project  we  continuously  investigated  further  application  possibilities  of 
the  newly  developed  methods.  Some  application  possibilities  were  already  mentioned, 
another  new  application  was  also  investigated  namely  the  application  of  a  the  generic 
optimal  controller  structure  for  raw  material  blending  [26]. 


General  remarks 

An  important  recognition  of  our  group  and  its  research  activity  is  the  invitation  for 
the  principal  investigator.  Prof.  L.  Keviczky  to  deliver  a  plenary  paper  on  our  recent 
results  at  the  IF  AC  Symposium  on  "Adaptive  Systems  in  Control  and  Signal  Processing" 
to  be  held  in  Budapest,  June  14-16,  1995. 

It  is  definitely  worth  mentioning  that  Dr.  Gary  Balas  spent  two  weeks  in  our 
Laboratory  last  summer  as  a  visiting  professor  from  the  The  Center  for  Control 
Sciences  and  Dynamical  Systems,  University  of  Minnesota,  which  is  our  partner 
from  the  US.  (The  first  version  of  our  application  was  a  joint  proposal,  but  later  it  was 
splited  to  a  Hungarian  and  a  US  one.)  His  visit  played  an  important  role  to  coordinate 
the  research  efforts  done  separately  and  to  present  jointly  the  scientific  achievements. 
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APPENDIX:  SELECTED  SCIENTIFIC  WORKS 


COMBINED  IDENTIFICATION  AND  CONTROL:  ANOTHER  WAY 
HOW  TO  INCREASE  ROBUSTNESS  OF  PID  REGULATORS 
CRITERIA  IN  FREQUENCY  DOMAIN 


APPLICATION  OF  A  NEW  GENERIC  OPTIMAL  CONTROLLER 
STRUCTURE  FOR  RAW  MATERIAL  BLENDING 

(See  in  separate  volume  ! ! !) 
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